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Abstract 

In  the  present  study,  we  propose  a  novel  configuration  of  partially  blocked  fuel  channel  with  baffle  plates  transversely  inserted  in  the 
channel.  The  effects  of  the  blockage  with  various  gap  ratios  and  numbers  of  the  baffle  plates,  the  fuel  flow  Reynolds  number  and  the  porosity 
of  the  diffusion  layer  on  the  reactant  gas  transport,  and  the  pressure  drop  across  the  channel  length  are  explored.  A  two-dimensional  model 
of  one-side  porous-walled  channel  with  the  presence  of  baffle  plates  is  considered.  An  approximate  analysis  is  performed  first  to  examine 
the  qualitative  natures  of  the  flow  in  the  gap  region  and  then  numerical  simulations  for  the  parametric  study  of  the  reactant  gas  transport 
in  a  half-cell  model  are  carried  out.  The  velocity  field,  the  oxygen  mass  flux/flow  rate  reaching  the  catalyst  layer,  the  concentration  of  the 
water  vapor  produced,  the  fraction  of  the  fuel  gas  entering  the  diffusion  layer,  and  the  pressure  drop  at  various  conditions  are  analyzed.  The 
results  reveal  that  reducing  the  gap  size  and/or  increasing  the  baffle  number  to  enhance  the  reactant  gas  transport  results  in  a  penalty  of  high 
pressure-loss.  Among  the  parameters  considered  in  the  present  work,  the  gap  ratio  has  the  most  remarkable  impact  on  the  variation  of  the 
pressure  drop.  Very  high  pressure-loss  can  be  generated  due  to  high  flow  resistance  at  a  low  gap  ratio  combined  with  a  low  gas  diffusion  layer 
(GDL)  porosity.  With  the  consideration  of  both  high  performance  and  reasonable  pressure  drop,  the  present  results  disclose  that  designs  with 
the  baffle  gap  ratio  no  smaller  than  0.1,  number  of  baffle  plates  N=  3-5,  and  the  GDL  porosity  around  0.7  seem  quite  appropriate. 
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1.  Introduction 

Design  of  the  flow  channel  in  bipolar  plates  is  one  of  the 
key  factors  in  the  system  performance  issue  of  proton  ex¬ 
change  membrane  (PEM)  fuel  cell  systems.  The  section  ge¬ 
ometry  and  the  pattern  of  the  fuel  flow  channel  influence  the 
reactant  gas  transport,  water  management,  and  the  efficiency 
of  the  utilization  of  the  fuel.  If  the  behaviors  of  the  exhaust 
gas  and  waste  water  are  out  of  control,  the  supply  of  oxy¬ 
gen  will  be  prohibited,  which  in  turn  may  cause  overheat  and 
dryout  of  the  fuel  cell  system.  This  unfavorable  situation  will 
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lead  to  the  increase  in  the  inter-resistance  and  deterioration 
in  the  system  performance.  On  the  contrary,  an  appropriate 
design  of  fuel  flow  channel  is  beneficial  to  the  reactant  gas 
transport  as  well  as  the  water  management. 

In  the  past  years,  therefore,  lots  of  efforts  have  been  en¬ 
deavored  to  novel  flow  channels  in  PEM  fuel  cell  systems 
such  as  the  arrangements  of  serpentine  channel,  multiple 
channels  in  parallel  type,  interdigited  channels,  etc.  Since  the 
fuel  channel  has  gas  diffusion  layer  (GDL)  as  a  side-wall,  the 
GDL  morphology  may  affect  the  reactant  gas  transport  from 
the  channel  to  the  catalyst  surface  and  is  usually  included 
in  dealing  with  the  gas  transport.  The  geometric  parameters 
have  been  analyzed  in  order  to  optimize  the  fuel  cell  system 
performance.  The  effects  of  the  rib  spacing  of  the  fuel  chan- 
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Nomenclature 

a  coefficient 

A,  B,  C  constants 
Cf  quadratic  drag  factor 

D\  mass  diffusivity  of  oxygen  (m2  s-1) 

£>2  mass  diffusivity  of  vapor  (m2  s-1) 

£a  activation  energy  ( J  mol- 1 ) 

h\  flow  channel  height  (m) 

h2  gap  between  baffle  and  diffusion  layer  (m) 

/?3  thickness  of  diffusion  layer  (m) 

k  permeability  (m2) 

L  flow  channel  length  (m) 

L2  baffle  plate  thickness  (m) 

m  mass  flow  rate  (kg  s-1) 

fn'[  molar  mass  flux  of  oxygen  (kg  m-2  s-1) 

M  molar  mass  (g  mol- 1 ) 

N  number  of  baffle  plates 

P  pressure  (N  m-2) 

R  gas  constant  ( J  mol- 1  K- 1 ) 

Re  Reynolds  number,  Uoh\/v 

t  time  (s) 

T  temperature  (K) 

U  velocity  in  v-direction  (ms-1) 

Uq  average  velocity  at  inlet  (ms-1) 

U  average  velocity  at  center  of  the  baffle  gap, 
yo*i/(*2  +  *3)(ms-1) 

V  velocity  in v-direclion  (ms-1) 

W\  molar  mass  fraction  of  oxygen 

W2  molar  mass  fraction  of  vapor 

X,  Y  cartesian  coordinates  (m) 

Greek  letters 

a  Beavers-Joseph  constant 

s  porosity 

X  gap  ratio,  fe/Zq 

v  kinematic  viscosity  (m2  s-1) 

p  density  (kg  m- 3 ) 

Subscripts 
0  inlet 

1  oxygen 

2  vapor 

b  boundary 

c  channel 

d  porous  diffusion  layer 

eff  effective 

i  ith  species 


nels  and  the  GDL  thickness  on  the  cell  performance  were 
studied  by  West  and  Fuller  [1].  By  conducting  experiments 
with  a  model  of  hydrophobilicity,  Jordan  et  al.  [2,3]  stud¬ 
ied  and  explained  the  effects  of  the  GDL  morphology  on  the 


cell  performance.  Gurau  et  al.  [4]  proposed  a  mathematical 
model  to  obtain  an  analytic  solution  of  the  mass  transport  of 
reactant  gas  in  a  half-cell,  in  which  the  effects  of  the  porosity 
and  the  tortuosity  of  the  gas  diffusion  layer  and  the  catalyst 
layer  were  explored. 

Recently,  numerical  computations  for  details  of  the  flow 
distribution  and  fuel  gas  diffusion  have  been  performed.  In 
their  simulation,  Hontanon  et  al.  [5]  focused  their  attention  on 
the  flow  distributor  to  search  for  optimization  of  the  fuel  flow 
characteristics.  They  employed  the  Navier-Stokes  equations 
and  the  Michaelis-Menten  type  two-step  kinetics  model  to 
study  the  performance  of  the  grooved  plate  and  porous  flow 
distributors  at  anode.  The  predictions  revealed  that  the  fuel 
consumption  increases  with  decreasing  permeability  of  the 
flow  distributor.  With  the  consideration  of  the  reactant  gas 
utilization  under  the  same  pressure  drop,  they  also  found  that 
the  porous  materials  are  more  advantageous  than  the  grooved 
plates.  Dutta  et  al.  [6]  performed  three-dimensional  compu¬ 
tations  for  fuel  transport  in  fuel  cells  of  serpentine  channels 
with  the  electrochemistry  modeled  by  the  mass  source/sink 
terms.  Kee  et  al.  [7]  proposed  a  generalized  computational 
model  for  the  mass  and  momentum  transport  in  channel  net¬ 
works  of  typical  planar  fuel  cell  layers/stacks.  Numerical  sim¬ 
ulation  was  performed  for  a  PEM  fuel  cell  with  a  metal  foam 
in  the  flow  field  of  the  bipolar  plates  [8]  and  the  results  demon¬ 
strated  that  the  use  of  metal  foam  renders  the  local  current 
density  distribution  more  uniform.  Yan  et  al.  [9]  conducted 
a  numerical  study  of  the  effects  of  fuel  channel  width  and 
GDL  porosity  on  the  cross-cell  transport  of  reactant  gas  and 
the  performance  of  a  PEM  fuel  cell.  Most  recently,  GDL  char¬ 
acteristics  with  various  thickness  and  porosity  fabricated  by 
using  different  manufacturing  methods  were  studied  by  nu¬ 
merical  simulation  and  experimental  measurements  [10]. 

From  the  above  brief  survey,  it  is  found  that  the  fuel  flow 
channel  design  is  very  significant  to  the  performance  of  the 
PEM  fuel  cells  and  various  configurations  have  been  stud¬ 
ied  previously.  The  interdigited  channel  with  fuel  flow  fully 
blocked  can  result  in  a  satisfactory  performance  in  reactant 
gas  transport.  In  the  previous  investigations,  however,  the 
penalty  of  high  pressure-loss  was  not  discussed  in  detail.  The 
objective  of  the  present  work  is  to  (1)  propose  a  novel  config¬ 
uration  of  partially  blocked  fuel  channel ;  (2)  explore  the  in¬ 
fluences  of  the  blockage  with  various  gap  ratios  and  numbers 
of  the  baffle  plates  on  the  reactant  gas  transport ;  (3)  evaluate 
the  pressure  drop  in  the  presence  of  various  blockage  effects ; 
and  (4)  examine  the  influences  of  fuel  flow  Reynolds  number 
and  porosity  of  the  gas  diffusion  layer.  A  two-dimensional 
model  of  one-side  porous-walled  channel  with  baffle  plate(s) 
in  it  is  considered.  The  hydrodynamic  nature  is  the  main 
theme  of  the  study  under  the  assumption  of  fast  reaction  at 
the  surface  of  the  catalyst  layer.  To  qualitatively  examine  the 
hydrodynamic  natures  and  used  as  a  limiting  case  for  nu¬ 
merical  verification,  an  approximate  analysis  of  the  flow  in 
gap  region  beneath  a  baffle  plate  is  performed  first.  Then  the 
numerical  simulation  is  carried  out  for  parametric  study  of 
the  mass/momentum  transport  in  the  whole  model  by  utiliz- 
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ing  the  commercial  software  CFDRC-ACE+.  The  velocity 
field,  the  oxygen  mass  flux/flow  rate  reaching  at  the  catalyst 
layer,  the  concentration  of  water  vapor  produced,  the  mass 
flow  fraction  of  the  fuel  gas  entering  the  diffusion  layer,  and 
the  pressure  drop  at  various  conditions  are  analyzed. 


2.  Half-cell  model  for  PEM  fuel  cells 


2.1.  Governing  equations 


Fig.  1  shows  the  present  two-dimensional  half-cell  model 
of  PEMFC  systems  for  the  fuel  gas  mass  and  momentum 
transport  in  the  flow  channel  and  the  porous  GDL.  The  fluid 
in  the  channel  is  air  at  inlet  pressure  Po  and  velocity  Uq.  The 
operating  temperature  is  set  to  be  To  =  333  K  in  the  simula¬ 
tion.  The  channel  length  is  L  =  L\  +L2  +  L3  with  L2  standing 
for  the  thickness  of  the  baffle  plate  and  L\  and  L3  for  the 
upstream  and  downstream  channel  lengths.  The  parameters 
h\,  Ii2  and  denote  flow  channel  height,  gap  size  and  GDL 
thickness,  respectively.  To  characterize  the  blockage  effect, 
a  dimensionless  parameter  named  gap  ratio,  X  =  h2/h\,  is  de¬ 
fined.  The  values  of  X  =  0  and  1  indicate  the  fully  blocked  and 
block- free  cases,  respectively,  and  the  values  in  between  mea¬ 
sure  various  levels  of  blockage.  The  flow  is  assumed  steady, 
laminar  and  incompressible,  the  porous  layer  is  of  isotropic 
and  uniform  material  with  the  porosity  e  and  the  permeabil¬ 
ity  k.  With  the  assumption  of  fast  reaction,  i.e.,  the  chemical 
reaction  time  scale  smaller  than  the  transport  time  scale,  the 
hydrodynamic  nature  plays  a  dominant  role  in  the  process 
and  the  chemical  reaction  is  expressed  simply  as  a  boundary 
condition  at  the  surface  of  the  catalyst  layer.  Since  the  cata¬ 
lyst  layer  is  extremely  thin  and,  therefore,  is  simply  regarded 
as  a  boundary  condition,  for  which  the  reactant  gas  (air  in 
this  case)  is  exhausted  for  perfect  reaction.  Under  the  above 
assumptions,  the  governing  equations  of  mass,  momentum 
and  species  balance  can  be  depicted  as  follows  [11,12]. 

(A)  In  the  flow  channel: 


du  dv 

dX  +  ~dY 


dU  dU 
U —  +  V  — 
8X  dY 


1  dP  (d2U  d2U\ 
pdX  +  V  V9^2  +  ~SY2) 


(1) 

(2) 


dV  dV 
U —  +  V —  = 
dX  dY 


1  3  P  Yd2V  d2V\ 
~pdY  +  V \3Z2  +  dY2) 


dWi  3  Wi 

U — -  +  V — -  =  Di 
dX  dY 


3  2Wi  d2Wi\ 
dX2  +  3  Y2  ) 


(3) 

(4) 


where  Wi  and  DL  denote  mass  fraction  and  diffusivity,  respec¬ 
tively,  the  subscript  i  =  1  stands  for  oxygen  and  i  =  2  for  water 
vapor. 

(B)  In  the  gas  diffusion  layer: 


du  dv 

dX  +  dY 


(5) 


'  dU  dU 
U —  +  V  — 
,  dX  dY 


s  dP  ( d2U  d2U\ 

~  ~~pJx  +  vs\dx2  +  Jy2) 

_v_^u_s^pUvu^V2  (6) 

k  v& 


<  dV  dv 

U —  +  V — 
,  dx  dY 


sdP  fd2v  d2v\ 
~~pdY+VS  +  W?) 

_vJy  (7) 

k 


m  dw,\  /d2Wi  d2Wi\ 

- h  V -  I  =  I  - —  4 - —  I 


/  3 Wi  dWj  \ 

U - -  +  V - - 

V  dX  3  Y  ) 


V  dX2  +  dY2  I 


(8) 


In  Eqs.  (6)  and  (7),  Cp  is  the  quadratic  drag  factor,  and 
Blake-Kozeny  correlation  [13]  is  adopted  for  the  relationship 
between  the  porosity  s  and  the  permeability  k ,  viz., 


Dd  = 


LI) 

\  £3  1 

V150/ 

L(i  —  e)2 J 

6V' 

% 

(9) 

(10) 


In  the  above  equations,  the  quantities  with  the  subscript  d 
denote  the  quantities  in  the  porous  diffusion  layer.  The  ge¬ 
ometry  parameter  V^/S'd  is  the  volume- to- surface  ratio  of  the 
GDL.  The  last  two  terms  in  Eqs.  (6)  and  (7)  are  the  drag  forces 
stemmed  from  the  presence  of  the  porous  wall  and  may  re¬ 
sult  in  an  increase  in  pressure  drop.  The  parameter  A;eff  is 


Fig.  1.  PEMFC  half-cell  model  of  partially  blocked  flow  channel. 
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the  effective  diffusivity  of  the  ith  gas  component  in  GDL  and 
the  Bruggemann  correlation  [14], 

A'.eff  =  Dj  ■  Sl  5  (11) 

is  adopted. 

2.2.  Boundary  conditions 

As  that  shown  in  Fig.  1,  the  conditions  at  the  boundaries 
designated  by  B1-B12  are  listed  as  follows: 

At  Bl: 

U  =  U0,  V  =  0,  P  =  Pq, 

Wi  =  Wifi,  W2  =  W2,o  (12a) 


U  (m/s) 


At  B2,  B6,  B8,  B12: 


dWi 

U=V=  — -  =  0 

dx 

At  B3,  B5,  B7: 
dWi 

U=V=  — -  = 0 

ay 

At  B4: 

U=V  =  Wi=  0, 
At  B9: 


dU 

~dY 

dWi 

~W 


=  £  - 


Y=ht 


du 

ay 


Y=h~ 


Y=h 


dWi 
+  ~S~dY 


Y=h~ 


(12b) 


(12c) 

(12d) 


dV 

~dY 


=  £  ■ 


Y=h3 


dV 

ay 


Y=h~ 


(12e) 


Fig.  2.  Grid-dependence  test  by  examining  the  axial  velocity  distributions 
at*/L  =  0.25. 


test  was  performed  to  validate  the  computations  and  exam¬ 
ine  the  grid-dependence  of  the  results.  By  employing  a  base¬ 
line  case  of  Reynolds  number  Re-  50,  one  baffle  plate  N  =  1, 
gap  ratio  X  =  0.05,  and  porosity  £  =  0.4,  computations  on  the 
four  grid  systems  with  number  of  grid  points  (length  x  width) 
215  x  30, 430  x  60, 645  x  60,  and  645  x  90  were  performed. 
The  axial  velocity  profiles  at  X/L  =  0.25  are  plotted  and  com¬ 
pared  in  Fig.  2.  It  can  be  observed  that  the  grid  of  430  x  60 
points  is  appropriate  with  the  considerations  of  accuracy  and 
time-saving. 


3.  Approximate  solution  of  flow  in  a  narrow  gap 
region  with  one-side  porous  wall 


At  B10: 

U=V  =  Wi=  0, 

d_Yll  =  oB'e~(E*'RT)  !  RT  x  x  ^5 

dY  P  y  2txM\  £>2,eff  9 Y 

At  B 1 1 : 


(120 


Fig.  1  shows  the  two-dimensional  flow  channel  model  with 
the  narrow  gap  region  —h 3  <  Y  <  h2  beneath  the  baffle  plate. 
To  facilitate  an  approximate  analysis,  assume  that  Darcy  law 
[16]  can  be  used  to  describe  the  fluid  flow  in  the  porous  GDL. 
The  governing  equation  of  the  flow  in  this  composite  channel 
configuration  can  be  written  as 

(A)  In  the  gap  region  of  the  flow  channel  0  <  Y <  h2\ 


dU  _  dV  _  dWi 
dx  ~  dx  ~  ~dx 


(12g) 


Since  the  fast  and  complete  reaction  of  oxygen  is  assumed 
at  the  interface  of  the  GDL  and  the  catalyst  layer,  the  relation 
of  the  oxygen  and  water  vapor  can  be  described  by  Eq.  (12f), 
where  B'  is  a  constant,  Ea  is  the  activation  energy,  R  is  the 
gas  constant,  and  M\  is  the  molar  weight  of  oxygen. 


V  •  V  =  0; 

dv  1  9 

—  +  V  •  XV  = - VF  +  vV2F. 

dt  p 

(B)  In  the  diffusion  layer  —  h2  <  Y  <  h2: 

V  •  Vd  =  0; 


(13) 

(14) 


(15) 


2.3.  Numerical  procedure 

The  system  of  the  governing  equations  described  above  is 
solved  utilizing  the  commercial  code  CFD-ACE+  based  on 
the  SIMPLE  Algorithm  [15]  for  convection-diffusion  prob¬ 
lems.  Before  mass  data  generation,  a  preliminary  numerical 


IdVd 

£  dt 


-lypd-V-V d, 

p  k 


The  related  boundary  conditions  are: 


At  Y  =  h2  :  V  =  (U,  V)  =  0. 


(16) 


(17) 
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dU  £ 

At  y  =  o :  —  =  —(u  -  udy,  v  =  yd; 

dY 

dv 

Pd  —  P  - h  2 [Ji —  —  0. 

^ dY 


(18) 


At  Y  =  -h3  :  Vd  =  0. 


(19) 


For  very  narrow  gap,  i.e.,  X  1,  the  flow  can  be  assumed 
fully-developed  in  the  narrow  gap.  The  above  momentum 
equation  can  be  reduced  to  the  form  of 


d2U  _  1  dP  _ 

_Vd YI~~~pdX~ 


(20) 


Integrating  the  last  equation  twice  gives  the  following  so¬ 
lution: 


AY2 

U  =  -—  +  BY 
2v 


C. 


(21) 


From  Eq.  (16)  with  —  (1/ p)dP/dX  =  A,  the  fluid  velocity 
in  GDL  becomes 


Ud 


-A. 

V 


(22) 


With  the  above  boundary  conditions,  the  constants  B  and 
C  in  Eq.  (21)  can  be  determined, 


Ah2 

B  =  — -L, 
v 

in  which 


C  = 


Ahi 


-M, 


(23) 


L  = 


ad 


M  = 


2(5  T  ad) 
5 


+ 


+ 


a25 
5  T  ad 
a82 


a8 

T’ 


2(5  +  ad)  d(8  +  ad) 


„  -  A 

^  i  ’ 

h3 


v  h2 
J  =  —. 
h2 

(24) 


The  fluid  velocity  in  the  channel  can  be  expressed  as 
hi  A  (  Y2  \ 

U(Y)=  1-—  +  LY  +  h2M)  .  (25) 

The  fraction  of  the  mass  flow  rate  through  the  diffusion 
layer,  mj,  with  respect  to  the  total  mass  flow  rate,  m0 ,  is  also 
an  important  parameter  in  design,  and  can  be  evaluated  by 
combining  Eqs.  (21)  and  (24),  viz. 

Wd  _  md  _ _ kh3 _  (2(f) 

m0~  md  +  mc~  kh3  +  h\[-  (1/6)  +  (L/2)  +  M] 

The  parameter  mc  stands  for  the  mass  flow  rate  through  the 
flow  channel.  With  the  relation  of  h2=kh\,  Eq.  (26)  can  be 
written  as 

md  md  kh2 

m0  ~  md  +  mc  ~  kh3  +  X3h\[-  (1/6)  +  (L/2)  +  M] 

(27) 


4.  Results  and  discussion 

4.1.  Effects  of  GDL  porosity  onflow  and  mass  transfer 

In  the  following  discussion,  a  baseline  case  of  Re  =  50, 
X  =  0.05,  and  e  =  0.4  is  defined  in  the  parametric  analysis.  To 
examine  effects  of  the  porosity,  the  values  of  s  =  0. 1,  0.4, 0.7, 
and  0.9  are  considered.  For  the  typical  cases  of  Re  =  50,  N=  1 , 
and  X  =  0.05  in  Fig.  3,  the  close-up  views  of  the  velocity  vector 
fields  in  the  baffle-gap-GDL  region  of  the  present  composite 
flow  passage  clearly  demonstrate  the  increase  in  mass  flow 
rate  through  GDL  of  high  porosity.  At  a  low  porosity,  the  gas 
flow  moves  downstream  mainly  through  the  gap  between  the 
baffle  and  the  GDL  at  a  high  speed,  thereby  a  recirculation 
appears  behind  the  baffle.  This  recirculation  bubble  may  be 
reduced  in  size  by  increasing  the  GDL  porosity.  For  example, 
with  the  flows  in  the  cases  of  s  =  0.4  and  0.7  compared  with 
that  of  £  =  0.1,  the  lower  velocity  at  the  exit  of  the  narrow 
gap  and  the  pushing  effect  of  the  flow  out  of  the  GDL  both 
tend  to  deflect  the  flow  behind  the  block  upward  and  reduce 
the  size  of  the  recirculation  bubble.  It  is  most  remarkable  in 
the  case  of  £  =  0.7.  However,  as  the  porosity  approaches  a 
value  close  to  unity,  e.g.,  £  =  0.9,  the  flow  velocity  slip  at  the 
interface,  y  =  0,  can  be  alleviated.  The  deflecting  or  pushing 
effect  stemmed  from  the  gas  flow  out  of  the  GDL  becomes 
weak  and  the  size  of  the  recirculation  bubble  can  be  enlarged 
to  a  certain  degree.  These  flow  characteristics  may  be  related 
to  the  pressure  drop  of  the  channel  flow. 


(a) 


Fig.  3.  Close-up  views  of  the  velocity  vector  field  near  the  gap  of  the  baffle 
plate  in  the  baffle-gap-GDL  region  of  the  present  composite  flow  passage  of 
one  baffle  at  Re  =  50  and  X  =  0.05. 
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Y+h3 

h2+h3 


(c) 


Fig.  4.  Comparisons  of  numerical  and  analytical  results  of  axial  velocity  distributions  in  the  gap  region  and  porous  GDL  layer  at  Re  =  50  and  A.  =  0.05.  (a) 
e  =  0. 1 ;  (b)  e  =  0.4;  (c)  e  =  0.7;  (d)  £  =  0.9. 


At  Re  =  50  and  A  =  0.05,  Fig.  4  presents  comparisons  of 
numerical  and  analytical  results  of  axial  velocity  profiles  at 
the  midway  of  the  gap  region,  i.e.,  X  =  L\  +L2/2.  The  thin 
line  of  Y  ~0.85  indicates  the  interface  of  the  flow  chan¬ 
nel  and  the  GDL.  For  the  small  values  of  the  porosity,  e.g., 
s  =  0. 1 ,  most  of  the  fuel  gas  passes  axially  through  the  gap  re¬ 
gion.  However,  with  increasing  porosity,  more  gas  enters  the 
porous  diffusion  layer  for  the  lower  flow  resistance.  There¬ 
fore,  the  fluid  velocity  and  the  corresponding  mass  flow  rate 
in  the  GDL  increase.  Although  the  quantitative  deviations 
exist,  both  the  approximate  analysis  and  the  computational 
results  show  the  same  trend  and  the  agreement  at  low  porosity 
is  better.  The  Darcy  law  gives  a  uniform  flow  solution  the  in 
porous  layer,  which  is  a  good  approximation  only  for  the  case 
of  low  porosity.  The  analytical  approximation  is  restricted  by 
the  validation  of  the  Darcy  law  with  noticeable  discrepancy 
at  higher  porosity. 

Fig.  5  reveals  the  effects  of  GDL  porosity  on  the  concen¬ 
tration  contours  of  the  oxygen  and  water  vapor.  The  model 


under  consideration  is  the  half-cell  on  the  cathode  side,  where 
the  oxygen  and  the  water  vapor  concentrations,  as  shown  in 
Fig.  5(a)  and  (b),  are  the  major  concerns  in  the  mass  transport 
aspect.  Fig.  5(a)  shows  the  oxygen  contours  corresponding 
to  the  cases  in  Fig.  3.  In  the  case  of  the  low  porosity  s  =  0.1, 
most  of  the  oxygen  appears  in  the  flow  channel,  only  very 
little  oxygen  can  transport  through  the  GDL  and  reach  the 
catalyst  layer.  This  leads  to  a  low  reaction  rate  and  there¬ 
fore  a  low  system  performance.  As  the  porosity  increases, 
e.g.  s  =  0.4,  0.7  and  0.9,  more  oxygen  enters  the  porous  layer 
and  the  oxygen  concentration  decreases  obviously  along  the 
channel  length.  The  reaction  at  the  catalyst  surface  is  en¬ 
hanced  with  more  fuel  supplied.  In  the  region  beneath  the 
baffle  plate,  a  large  amount  of  oxygen  is  forced  to  turn  into 
the  porous  GDL  and  strong  reaction  and  higher  consumption 
of  oxygen  occur  in  this  region. 

Water  is  the  product  of  the  reaction  at  the  surface  of  the 
catalyst  layer,  therefore,  the  maximum  water  vapor  concen¬ 
tration  can  be  observed  at  the  interface  of  the  GDL  and  the 
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(a)  X(mm)  (b)  X(mm) 


Fig.  5.  Effects  of  porosity  on  the  local  oxygen  and  water  vapor  concentration  in  the  baffle-gap-GDL  region  of  the  present  composite  flow  passage  of  one-baffle 
at  Re  =  50  and  X  =  0.05.  (a)  Oxygen;  (b)  water  vapor. 


catalyst  layer  where  the  oxygen  is  consumed.  For  low  GDL 
porosity,  the  reactant  gas  transport  through  the  GDL  is  diffi¬ 
cult  and  the  reaction  is  weak.  However,  the  situation  can  be 
improved  by  increasing  the  porosity  of  the  GDL.  Fig.  5(b) 
demonstrates  an  increasing  water  vapor  transport  from  the 
catalyst  layer  to  the  flow  channel  with  porosity.  Therefore, 
high  porosity  is  beneficial  to  remove  the  excess  water.  In 
Fig.  5(a)  and  (b),  it  is  also  observed  that  the  concentrations 
of  the  oxygen  (reactant)  and  the  water  vapor  (product)  are  of 
an  approximately  complementary  relation.  The  higher  wa¬ 
ter  concentration  appears  at  the  location  where  there  is  less 
oxygen,  and  vice  versa. 

At  the  catalyst  surface,  the  effects  of  GDL  porosity  on 
the  axial  distribution  of  the  oxygen  mass  flux  are  displayed 
in  Fig.  6.  In  general  the  oxygen  decreases  along  the  ax¬ 
ial  length  but  with  an  abrupt  jump  in  the  region  beneath 
the  baffle  plate,  where  the  strong  forced  convection  en¬ 
hances  the  transport  of  the  oxygen  as  mentioned  previously. 
The  peak  of  the  mass  flux  jump  is  raised  up  in  the  cases 
of  high  porosity,  e.g.,  £  =  0.7  and  0.9.  Behind  the  baffle 
plate,  however,  the  oxygen  fluxes  for  high  porosity  become 
less  due  to  the  stronger  reaction  with  fuel  depletion  occurs 
upstream. 


Fig.  6.  Effects  of  porosity  on  the  axial  distribution  of  oxygen  mass  fluxes  at 
the  surface  of  the  catalyst  layer. 
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Fig.  7.  Axial  velocity  distributions  in  the  gap  region  and  porous  GDL  layer  at 
and  s  =  0.4;  (b)  effects  of  Reynolds  number  at  k  =  0.05  and  s  =  0.4. 

4.2.  Effects  of  gap  size  and  Reynolds  number 

It  can  be  expected  that  the  change  in  flow  area  of  the  fuel 
channel  may  have  a  considerable  influence  on  the  mass  trans¬ 
port  in  the  flow  channel  and  the  GDL  as  well.  The  blocked 
channels  with  baffle  of  gap  ratios  A  =  0, 0.005, 0.05, 0. 1  and  1 
are  investigated  in  the  present  work.  For  the  cases  of  Re  =  50 
and  £  =  0.4,  Fig.  7(a)  shows  the  blockage  effects  on  the  veloc¬ 
ity  profiles  at  the  midway  location  X  =  L\  +  Lff2  of  the  gap 
region  between  the  baffle  and  the  GDL.  An  abrupt  change 
of  velocity  occurs  at  the  interface  except  in  the  case  of  fully 
blocked  channel  (A.  =  0),  in  which  all  of  the  fluid  is  forced  to 
enter  and  flow  through  the  GDL.  As  the  flow  channel  is  par¬ 
tially  blocked  (0  <  k  <  1),  however,  part  of  the  fluid  is  allowed 
to  move  downstream  through  the  gap  and  the  rest  is  forced 
into  the  GDL.  With  the  increase  in  the  gap  height,  the  mass 
flow  rate  through  the  gap  increase  but  the  gap  flow  velocity 
may  reduce  with  the  enlarged  flow  area.  Due  to  the  counter¬ 
action  of  these  two  effects,  the  gap  flow  velocity  variation 
with  the  gap  ration  k  is  not  monotonic.  The  data  shown  in 
Fig.  7  demonstrates  that  the  axial  velocity  at  k  =  0.05  has  a 
maximum  among  the  cases  considered.  While  the  Reynolds 
number  has  only  minor  influence  on  the  axial  velocity  profile 
in  the  gap  region  as  shown  in  Fig.  7(b). 

The  amount  of  reactant  gas  entering  the  porous  GDL  is 
one  of  the  most  interested  and  significant  flow  quantities. 
Fig.  8  presents  the  ratio  of  the  mass  flow  rate  in  porous  GDL 
to  the  total  mass  flow  rate,  mp0rousMtotai>  at  various  flow  and 
geometric  conditions.  It  is  observed  that,  in  general,  the  nu¬ 
merical  predictions  of  rh porous  M  total  are  of  reasonable  agree¬ 
ment  with  the  analytic  approximations.  From  the  results,  it 
is  also  found  that  this  fraction  of  the  mass  flow  rate  in  the 
GDL  reduces  as  the  gap  ratio  increases.  For  a  fixed  value  of 
A,  the  fraction  of  the  reactant  gas  transport  in  the  GDL  in¬ 
creases  with  increasing  Re  since  a  higher  pressure  gradient, 
and  thus  the  higher  driving  force,  can  be  built  in  the  presence 


the  central  section  of  the  blocked  gap  region,  (a)  Effects  of  gap  size  at  Re  =  50 


of  the  baffle  plate.  The  Re  effects  diminish  at  very  small  gap 
ratio  (A  — >►  0)  and  very  large  ones  (A  — ►  1)  as  well.  Among 
the  cases  studied,  the  Re  effect  at  A  =  0.2  is  most  noticeable. 

Fig.  9  discloses  effects  of  gap  ratio  and  Reynolds  number 
on  the  axial  distribution  of  the  oxygen  mass  flux,  m",  on 
the  catalyst  surface.  In  Fig.  9(a),  the  oxygen  mass  flux  m![ 
is  shown  at  Re  =  50  and  s  =  0.4.  The  distribution  of  fn'[  is 
characterized  by  a  descending  trend  along  the  channel  length 
but  has  a  peak  value  at  the  location  beneath  the  baffle  plate. 
Larger  amount  of  reactant  mass  is  forced  towards  the  catalyst 
surface  by  the  presence  of  the  baffle  plate,  and  the  smaller 
the  gap  size,  the  higher  the  mass  flux  is.  However,  changes 
of  m![  appearing  in  the  upstream  and  downstream  regions 
of  the  baffle  plate  are  moderate.  The  oxygen  mass  flux  in 
the  channel  can  be  enhanced  by  increasing  Re  which  gives  a 
strong  mass  convection  in  the  channel  main  flow,  as  shown 
in  Fig.  9(b). 


Fig.  8.  Normalized  mass  flow  rates  in  porous  layer  at  the  axial  location  of 
gap  region  center  and  at  various  conditions. 
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Fig.  9.  Oxygen  mass  flux  distributions  at  the  surface  of  the  catalyst  layer,  (a)  Effects  of  gap  size  at  Re  =  50  and  £  =  0.4;  (b)  effects  of  Reynolds  number  at 
.A.  =  0.05  and  £  =  0.4. 


4.3.  Effects  of  number  of  baffle  plates  on  mass  transport 

Number  of  the  baffle  plate  is  one  of  the  significant 
parameters  in  the  present  analysis  of  transport  phenomena. 
Fig.  10  shows  the  distributions  of  the  oxygen  flux  mf[  along 
the  interface  between  the  catalyst  layer  and  the  GDL  at 
Re-  50,  £  =  0.7,  and  A  =  0,  0.05,  and  0.1.  Locally,  the  mass 
flux  is  of  peak  values  at  the  locations  corresponding  to 
the  baffle  plates  and  the  peak  values  decay  downstream. 
With  various  numbers  of  the  baffle  plates,  the  cases  with 
small  gap  ratio  A  generally  have  higher  peak  values  of  the 
oxygen  mass  flux  under  the  first  baffle.  The  total  mass  flow 
rate  of  the  oxygen  reaching  the  catalyst  layer,  m  i,  can  be 


evaluated  by  the  integration  rii\  =  f  \m"  dx,  where  x  =  X/L. 
Fig.  11  presents  the  values  of  m\  for  the  cases  of  Re  =  50 
and  £  =  0.7  with  various  gap  ratio  A  and  baffle  number 
N.  Reduction  in  gap  ratio  causes  an  increase  in  m\  due 
to  blockage  effects  of  the  baffle  plates.  Also,  the  data  in 
Fig.  1 1  demonstrate  that  adding  baffle  plates  enhances  the 
oxygen  mass  transport  toward  the  catalyst  layer  through 
GDL. 

From  the  data  shown  in  Fig.  10,  it  is  obvious  that  the 
axial  distributions  of  the  oxygen  mass  flux  mf[  at  the  inter¬ 
face  of  the  GDL  and  the  catalyst  have  a  similar  pattern.  In 
Fig.  12,  we  plot  the  normalized  mass  flow  rate  mp0roUs/^totai 
at  the  midway  of  the  channel,  X/L  =  0.5,  under  conditions  of 


Fig.  10.  Axial  distributions  of  oxygen  mass  fluxes  at  Re  =  50,  £  =  0.7  and  various  values  of  gap  size  and  baffle  number,  (a)  N=  3;  (b)  N=  5;  (c)  N= 7;  (d)  N=  9. 
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Fig.  11.  Oxygen  mass  flow  rate  reaching  the  interface  of  GDL  and  catalyst 
layer  for  the  cases  of  Re  =  50  ,s  =  0.7  and  various  values  of  gap  size  and  baffle 
number. 

various  combinations  of  porosity,  gap  ratio,  and  number  of 
baffle  plates  to  demonstrate  the  similarity.  As  the  trend  shown 
in  Fig.  8,  at  a  larger  gap  ratio,  the  fraction  mp0rous/^total  is 
lower  since  most  of  the  fluid  moves  downstream  through  the 
gap  between  the  baffle  plates  and  the  GDL,  and  a  smaller 
amount  of  fluid  is  forced  to  pass  through  the  porous  GDL. 
On  the  other  hand,  the  mass  flux  fraction  increases  with  in¬ 
creasing  porosity  for  the  smaller  flow  resistance  in  the  GDL. 
However,  this  special  value  of  the  mass  flux  fraction  at  the  lo¬ 
cation  XI L  =  0.5  is  almost  independent  of  the  number  of  baffle 
plates. 

To  examine  the  influence  of  baffle  plate  number  on  the  lo¬ 
cal  distribution  of  water  vapor  concentration,  contours  at  the 
surface  of  the  catalyst  layer  for  the  cases  of  Re  =  50,  e  =  0.4, 
k  =  0.05,  and  N=  3,  5,  7,  and  9  are  shown  in  Fig.  13.  Under 
the  assumption  of  fast  and  complete  reaction  of  the  oxygen 
reaching  the  surface  of  the  catalyst  layer,  the  amount  of  the 
water  generated  directly  depends  on  the  local  mass  flux  of 
the  reactant  gas,  i.e.,  oxygen.  It  is  observed  that  high  water 
vapor  concentration  appears  at  the  locations  under  the  baffle 
plates  where  there  is  a  high  oxygen  mass  flux  as  shown  in 
Fig.  10. 


e 

Fig.  12.  Effects  of  baffle  plate  number,  gap  ratio  X  and  porosity  s  on  fraction 
of  the  mass  flow  rate  entering  the  GDL  at  X/L  =  0.5. 
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Fig.  13.  Effects  of  baffle  number  on  the  local  water  vapor  concentration  in 
the  baffle-gap-GDL  region  of  the  present  composite  flow  passage  at  Re  =  50, 
A.  =  0.05  and  £  =  0.4. 

4.4.  Pressure  drop  across  the  channel  length 

The  presence  of  the  baffle  plates  placed  normal  to  the  main 
flow  direction  partially  blocks  the  fuel  channel  and  forces 
most  of  the  fuel  to  go  through  the  porous  GDL  and  reach 
the  catalyst  interface.  This  arrangement  improves  the  fuel 
transport  rate  and,  therefore,  the  reaction  at  the  catalyst  layer 
can  be  enhanced.  However,  the  intrusion  of  the  baffle  plates 
causes  larger  pressure-loss  and  needs  higher  pumping  power 
for  delivery  of  the  fuel.  Fig.  14  presents  the  influences  of  baf¬ 
fle  gap  ratio,  main  flow  Reynolds  number,  and  the  number  of 
baffle  plates  on  the  dimensionless  pressure  drop  across  the 
fuel  channel,  A  F*  =  (Piniet  —  Pexit)/(pUQ/2).  From  the  data 
in  Fig.  14(a),  it  can  be  realized  that  the  pressure  drop  increases 
with  reducing  gap  size  and/or  porosity.  With  an  appropriate 
gap  between  the  baffle  and  the  GDL,  e.g.,  k  =  0.1  or  larger, 
the  pressure  drop  and  then  the  pumping  cost  needed  would 
be  considerably  reduced.  The  data  in  Fig.  14(b)  show  that 
the  dimensionless  pressure  drop  reduces  with  the  increasing 
main  flow  Reynolds  number  and  the  effects  in  the  cases  con¬ 
sidered  are  moderate.  This  trend  is  similar  to  the  theory  of  the 
laminar  channel  flow,  in  which  the  pressure  drop  is  inversely 
proportional  to  the  Reynolds  number.  Increasing  the  number 
of  the  baffle  plates  implies  the  stronger  blockage  effect  and 
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Fig.  14.  Pressure  drop  in  one-side  porous-walled  channel  at  various  con¬ 
ditions.  (a)  Effects  of  baffle  gap  size  for  Re- 50  and  N=  1;  (b)  effects  of 
Reynolds  number  for  N=  1  and  k  =  0.1;  and  (c)  effects  of  number  of  baffle 
plates  for  Re  =  50  and  X  =  0. 1 . 

higher  pressure-loss,  as  shown  in  Fig.  14(c).  From  the  present 
results,  it  is  found  that  extremely  high  pressure  loss  can  be 
generated  due  to  high  flow  resistance  at  a  low  gap  ratio  com¬ 
bined  with  a  low  GDL  porosity.  In  general,  the  pressure  drop 
in  the  fuel  channel  can  be  reduced  noticeably  with  increasing 
porosity  of  GDL  in  the  range  of  s  >  0.6. 


5.  Concluding  remarks 

In  the  present  study,  we  have  investigated  the  reactant  gas 
transport  in  a  PEM  fuel  cell  with  the  fuel  channel  partially 
blocked  by  installing  baffle  plates.  The  effects  of  the  baffle 
blockage  and  the  porosity  of  the  GDL  on  the  fuel  gas  flow 
and  mass  transfer  have  been  examined  in  detail.  Based  on  the 
present  results  and  analysis,  the  following  conclusions  can 
be  drawn. 


The  presence  of  the  baffle  plates  has  significant  impact  on 
the  flow  field  in  the  fuel  channel  and  the  gas  reactant  trans¬ 
port  through  the  GDL.  The  baffles  block  the  main  flow  in  the 
channel  and  force  more  fuel  gas  to  turn  to  the  GDL.  This  fact 
implies  an  enhancement  of  the  oxygen  mass  flux  and  the  water 
vapor  concentration  at  the  catalyst  surface,  especially  at  the 
positions  right  beneath  the  locations  of  the  baffle  plates.  The 
simulation  demonstrated  that  the  diffusion  layer  of  £  =  0.7 
is  beneficial  to  the  performance  of  the  reactant  gas  transport. 
Increasing  the  number  of  baffles  in  the  flow  channel,  the  oxy¬ 
gen  flux  in  the  catalyst  layer  under  the  baffle  would  become 
higher.  However,  reducing  gap  size  and/or  increasing  baf¬ 
fle  number  to  enhance  the  reactant  gas  transport  results  in 
a  penalty  of  high  pressure-loss.  Although  the  present  model 
without  three-dimensional  relieving  effect  may  exaggerate 
the  pressure-loss,  the  qualitative  trend  of  the  present  predic¬ 
tions  is  useful  for  evaluation  of  the  penalty  stemmed  from 
high  resistance  due  to  the  baffle  blockage  and  the  low  GDL 
porosity.  In  the  parameter  ranges  considered  in  the  present 
work,  the  gap  ratio  has  the  most  remarkable  impact  on  the 
variation  of  the  pressure  drop.  Extremely  high  pressure  loss 
can  be  generated  due  to  high  flow  resistance  at  a  low  gap  ratio 
combined  with  a  low  GDL  porosity  and  should  be  avoided. 
With  the  consideration  of  both  high  performance  and  reason¬ 
able  pressure  drop,  the  present  results  disclose  that  designs 
with  the  baffle  gap  ratio  no  smaller  than  0.1,  number  of  baffle 
plates  N=  3-5,  and  the  GDL  porosity  around  0.7  seem  quite 
appropriate. 
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